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Multi-dimensional core-collapse supernova simulations exhibit turbulence of large amplitude and 
over large scales. As neutrinos pass through the supernova mantle the turbulence is expected to 
modify their evolution compared to the case where the explosion is free of turbulence. In this paper 
we study this turbulence effect upon the neutrinos modelling the turbulence expected from multi- 
dimensional simulations by adding matter density fluctuations to density profiles taken from one- 
dimensional hydrodynamical simulations. We investigate the impact upon the supernova neutrino 
transition probabilities as a function of the neutrino mixing angle On and turbulence amplitude. In 
the high (H) resonant channel and with large ^13 values we find that turbulence is effectively two 
flavor for fluctuation amplitudes < 1% and have identified a new effect due to the combination of 
turbulence and multiple H resonances that leads to a sensitivity to fluctuations amplitudes as small 
as ~ 0.001%. At small values of beyond the range achievable in Earth based experiments, we 
find that turbulence leads to new flavor transient effects in the channel where the MSW H resonance 
occurs. Finally, we investigate large amplitude fluctuations which lead to three flavor effects due 
to broken HL factorization and significant non-resonant transitions and identify two non-resonant 
turbulence effects, one depending on the ^13, and the other independent of this angle and due to 
the low (L) MSW resonance. 
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I. INTRODUCTION 

Core-collapse supernova simulations have reached a 
high degree of sophistication and complexity in the last 
decade. A new picture is emerging that shows how 
several ingredients might be necessary to achieve suc- 
cessful explosions. These include fluid instabilities (the 
SASI mode), realistic nuclear networks, detailed neutrino 
transport, rotation, magnetic fields etc. [l|,[2|- The flavor 
evolution of the neutrino signal from the next Galactic 
supernova as a function of time and energy has the poten- 
tial to reveal invaluable information about the explosion 
of the star and is being thoroughly investigated for the 
case of spherically symmetric explosions. 

Fully understanding neutrino flavor conversion in the 
supernova environment and the complex interplay with 
the various features of the supernova dynamics is clearly 
a theoretical challenge. Nevertheless, there has been sig- 
nificant ongoing progress and recent developments have 
highlighted a variety of novel flavor conversion phenom- 
ena. This is due to the inclusion of both the neutrino- 
neutrino interaction and the use of dynamic density pro- 
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files taken from simulations. While the former induces 
collective effects the latter produces dynamic mul- 
tiple Mikheev-Smirnov-Wolfenstei n ( MSW) resonances 
which can lead to phase effects |9l-[l6l|. For a review see 
[13, [13 • The first calculation for the supernova neutrino 
signal implementing both of these effects together consis- 
tently is performed in [l^ while the impact of collective 
and dynamic MSW effects upon the diffuse supernova 
neutrino background can be found in [lO, HH . 

Multi-dimensional supernova simulations indicate 
strong turbulence, particularly in the region between the 
reverse and forward shocks, generated by non-radial fluid 
flow through the distorted shock surfaces [2^ - |2^ and 
other multi-dimensional phenomena. But due to the ab- 
sence of suitable multi-dimensional simulations impor- 
tant properties of the turbulence, such as its amplitude at 
late times, are unknown. Turbulence effects upon neutri- 
nos, not just in supernova, have been investigated in sev- 
eral works [25l - [28j . In [2^ neutrino oscillations in noisy 
media with ^-correlated fluctuations is considered with 
applications to a fluctuating matter density and magnetic 
flelds in the sun. In it is also found that depolariza- 
tion can also occur even if the fluctuations never cause 
the density to cross the neutrino resonance density i.e. 
they are off-resonance. The application to supernovae is 
studied in [soj with a similar model for the turbulence to 
determine the implications for neutrino flavor conversion 
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and supernova dynamics (r-process and shock reheating). 
In [3l| the authors use the same model of deha-correlated 
fluctuations for core-collapse supernovae, but now with 
profiles containing both front and reverse shocks. A sen- 
sitivity to the third neutrino mixing angle is shown in 
the channels where no MSW H resonance occurs (elec- 
tron neutrino channel in inverted hierarchy and electron 
anti- neutrino channel in normal hierarchy). Finally ar- 
guments were made in [sl] about the application of Kol- 
mogorov correlated fluctuations and linear profiles to su- 
pernova. Regardless of their differences, all these studies 
agree that the inclusion of matter density fluctuations in 
media has the possible consequence of producing com- 
pletely mixed states (a process called 'depolarization') if 
the amplitude is sufficiently large. At the same time all 
these investigations are based on the solutions of differen- 
tial equations for the averaged density matrix using just 
two neutrino flavors. The first full three flavor study of 
turbulence effects is performed in Ref.js^ where turbu- 
lence is superimposed upon profiles from one-dimensional 
hydrodynamical simulations using a Kolmogorov spec- 
trum for the density fluctuations. In contrast to the pre- 
vious work, here turbulence effects are studied by build- 
ing a statistical ensemble of instantiations for the neu- 
trino survival probability. It is also shown there that 
HL factorization might be broken in the limit of strong 
turbulence. But despite the acknowledgment that tur- 
bulence may be important the effect is not considered in 
the calculations for the neutrino signal. Thus it remains 
unclear if the presence of turbulence will remove the in- 
formation about the explosion imprinted in the neutrino 
signal and which is supposed to indicate the missing in- 
formation about neutrino mixing. 



The purpose of this paper is to revisit the subject of 
turbulence effects upon core-collapse supernova and ex- 
plore the effects of changing the angle 6*13 and the tur- 
bulence amplitude. As in Ref. [s^ we follow the neu- 
trino evolution in the star using matter density profiles 
from a one-dimensional hydrodynamical simulation with 
Kolmogorov correlated turbulence superimposed. Simi- 
larly we construct an ensemble of results computed using 
a phase-retaining 3-flavor integrator instead of determin- 
ing the average density matrix and survival probabilities. 
We then investigate the distributions of the probabilities 
and compute the mean and the variance in order to study 
how the effects we find evolve as a function of the param- 
eters. Our paper is structured as follows: in Section II 
the model and the numerical procedures are described. 
Section III presents the numerical results of turbulence 
effects on the survival probabilities first for the case of 
small amplitude turbulence and both large and small 613 
and then large amplitude fiuctuations which can cause 
three flavor mixing in both the H-resonant and non-H- 
rcsonant channel. Section IV contains our conclusions. 



II. THE MODEL 

Neutrinos evolve through matter according to the 
Schrocdinger-likc equation: 

^J^\^)=H\^P) = {K + V)\^/J). (I) 

In the absence of neutrino-neutrino interactions the 
Hamiltonian is composed of two terms: the vacuum term 
K, and the potential term V due to the matter. In the 
flavor basis the vacuum term K^^^^ depends upon two 
mass-square differences Smi2 and Sm^^ and on the pa- 
rameters of the Maki-Nakagawa-Sakata-Pontecorvo ma- 
trix which arc the three neutrino mixing angles usually 
labeled as 6*12, 9i3, 623 and three phases, one Dirac and 
two Majorana. For all results in this paper we set the os- 
cillation frequencies and angles to (Smjj = 8 x 10~^eV^, 
\6ml^\ = 3 X lO-^eV^ and by sin^ 26I12 = 0.83 and 
sin^ 2023 = 1 [IJ ■ We shall consider both normal and 
inverse hierarchies and various values for the unknown 
angle ^13 up to its hmit from Chooz [H, . We do not 
consider any effect coming from the Dirac CP phase S. 
For a discussion of its impact see [1, [s^l ■ 

The potential term, at some fixed time is diagonal in 
the flavor basis i.e. ^'^•^^(r) = diag{Ve{r),Vfj,{r),Vr{r)) 
with 14 (r) = \/2GFne{r) and rie(r) the electron den- 
sity. The potentials V^i(r) and Vr{r) are negligible com- 
pared to T4(r) and will be ignored hereafter. The tur- 
bulence we shall study in this paper enters into equation 
([T]) through rie so in order to compute the neutrino evo- 
lution one needs to provide a density profile from which 
we can construct Ve- Ideally a study of the effect of 
turbulence in supernova upon the neutrinos would use 
density profiles taken from successfull multi-dimensional 
core-collapse supernova simulations. While significant 
progress continues to be made in the sophistication of 
the simulations we are still far from being able to extract 
fiuctuation characteristics (scale, power spectra and am- 
plitude) from them at late times and in the outer re- 
gions of the supernova where our interest lies. There- 
fore one must adopt a pragmatic approach by using a 
'turbulence' free explosion and then adding the turbu- 
lence to it in some prescribed fashion [3l|, [SJ]- The scale, 
power spectrum and amplitude of the turbulence are thus 
variable. While this approach requires calibration before 
it can be applied to understanding turbulence in multi- 
dimensional supernova simulations, an exploration of the 
turbulence parameter space can discover general features 
of turbulence effects upon the neutrinos. There are also 
several advantages of this approach: first, the underlying 
profile is the same so that non-turbulent features - such 
as the shocks - are always located at the same positions, 
second the turbulence in the profiles of the members of 
an ensemble is uncorrelated, and third, if the turbulence 
has a zero expectation value then the average potential 
(V) is known exactly allowing us to compare the effects 
of turbulence to the situation when turbulence is absent. 

Now that we have outlined our approach we spell out 
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the details. We start with a profile taken from a spheri- 
cally symmetric simulation, i.e. a one dimensional model. 
We shall ensure that this one dimensional profile is also 
the average profile {V){r) once the turbulence is added. 
The one-dimensional density profile we employ is the 
4.5 s siiapshot taken from the hydrodynamical simulation 
of Ref.fS^l with Q = 3.36 x 10^^ erg. This particular simu- 
lation was chosen because it has the greatest resemblance 
with the 2D hydrodynamical simulations there and be- 
cause the density of the region between the shocks - into 
which we will add the turbulence - overlaps with the H 
resonance density for neutrino e nerg ies between 10 and 
80 MeV as shown in Calais et al. [2l[ . Although we show 
results for a given snapshot in time our findings are valid 
for the entire period the shock wave is in the H resonance 
region for neutrino energies of order 0(10 McV) i.e. from 
t ^ 2 s to t '-^ 10 s or so. Later the shocks will move to 
affect the L resonance but the neutrino flux will be so 
small at these times that turbulence effects may be well- 
nigh impossible to observe. The turbulence is included 
in the model by writing the potential term in Eq.([T]) as 

V{r) = {l + F{r)){V){r) (2) 

with F{r) a random field. The field F{r) is constructed 
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X ^ {An COS [kn (r - rr)] + Bn sin [kn (r - r^)]} 

n=l 

(3) 

for Tr < r < rs and is zero outside this range. In this 
equation the parameter sets the amplitude of the fluc- 
tuations. The two radii and are the positions of 
the reverse and forward shock respectively; the two tanh 
terms are included to suppress fluctuations close to the 
shocks and prevent discontinuities at r,, and r^, and the 
parameter A is a scale over which the fluctuations reach 
their extent size. We set A = 100 km. The second half 
of equation (O is the discrete Fourier representation of 
a random held. The members of the set of co-efficients 
{A} and {B} are independent standard Gaussian ran- 
dom variates with zero mean thus ensuring the vanishing 
expectation value of F. To determine the Nk fc's, A's and 
J5's for an instantiation of F the 'Randomization method' 
described in Ref.js^ is employed. The number of modes 
used throughout most of this paper is Nk = 100 but to 
assure the reader that this is sufficient we shall also com- 
pute one of our principal results using = 1000 so that 
the reader can observe that the results for Nk — 100 arc 
essentially similar. Finally, the wavenumber cutoff k^, is 
set to ki, = 7r/(rs — Vr) i.e. a wavelength twice the dis- 
tance between the shocks. The power spectrum E{k) we 
use is 



i?(fc) = (a-l)(^ 



(4) 



Throughout this work we shall adopt a Kolmogorov spec- 
trum where a = 5/3. Note that the power spectrum in 
Eq.(l4]) is very different from the one of the (5— correlated 
power spectra of Refs.fsoj and [3T| . 

Now that all terms in Eq.(IT]) are defined, we solve the 
three flavor neutrino evolution equation with the phase- 
retaining integrator of Ref. (40| . The quantities we are in- 
terested in calculating, are the neutrino and antincutrino 
survival and appearance probabilities, i.e. P{iyj — > Ui) 
and P{vj — > i^i), in the matter basis since the flux of 
neutrinos emerging from the star is given by these quan- 
tities, multiplied by the appropriate initial fluxes at the 
proto-neutron star. The relationship between the flavor 
and matter bases can be found in Kneller & McLaughlin 
[40j . Such probabilities can be formed from the wave- 
function for each initial pure state or from the elements 
of the 5-matrix. For brevity we shall often use the nota- 
tion Pij = Pii'j — > Vi) and Pij = P{vj — > Di). 



III. NUMERICAL RESULTS: TURBULENCE 
IMPACT ON THE PROBABILITIES 

Our aim is to explore the parameter space of our model 
in order to survey the various effects turbulence produces 
and derive the relationships between the survival and 
crossing probabilities as a function of the parameters. 
This represents a large undertaking, for now we shall 
limit our parameter space to the mixing angle ^13 and 
the fluctuation amplitude and shall reserve for a later 
date [4l| an exploration of the effects of changes in the 
turbulence power spectrum, the snapshot proflle and the 
turbulence location. In particular we shall not consider 
the effect of turbulence in front of the forward shock, 
an issue intimately tied to the observability of neutrino 
transformation effects in the neutrino burst from the next 
galactic supernova. 

We begin with the case of large 6'i3 and explore the ef- 
fects in the H resonant channel before considering smaller 
values of 6*13. We shall flnish with the effects of large am- 
plitude fluctuations that can break HL factorization and 
lead to large effects in the non-resonant channels. 



Effective two flavor mixing in the H resonant 
channel 



By following the neutrino as it progresses through the 
proflle for a given instantiation we discover that the effect 
of turbulence can span a range from the pcrturbativc to 
the strong. First, turbulence introduces many new H res- 
onances which lead to substantial changes from the evo- 
lution of the same neutrino through the proflle sans tur- 
bulence. We also observe the perturbative effect caused 
by non-resonant fluctuations i.e. those fluctuations of the 
density that approach, but do not intersect, the reso- 
nance density. These two limits, the strong and the per- 
turbative, can be seen in flgurc ^ where we show the 
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FIG. 1: One instantiation of the density profile (top panel) 
and tlie survival probability P22 = P(i'2 — ^ 1^2) for a E — 
25 MeV neutrino (bottom panel) as functions of the radius r. 
The turbulence region is identified by the vertical dashed lines 
in both panels, the H resonance density is the horizontal line 
in the upper panel determined using the two flavor formula. 
The underlying base profile is shown as the dashed line in 
the upper panel, the evolution through the base profile as the 
dashed line in the lower panel. The value of 613 is sin^ (26^13) = 
4 X 10"*, the hierarchy is normal, and the amplitude of the 
density fluctuations is = 0.1. 



evolution of the second matter eigcnstate survival prob- 
ability, P22 (lower panel), as a function of the radius r 
through one instantiation of the matter density profile 
(top panel). The results are for a 'large' value [42[ of ^13 
corresponding to sin^ 2^13 = 4x 10"'', a 25 MeV neutrino 
and the normal hierarchy. We see that vi^ithout fluctua- 
tions the neutrino experiences just three H resonances for 
this chosen profile and energy, and one is located outside 
the turbulence region. After the inclusion of turbulence 
the number of MSW resonances increases significantly. 
Since such resonances are semi-adiabatic for this value of 
013 the presence of the extra resonances due to turbu- 
lence causes the neutrino survival probability to undergo 
large transitions as the neutrino passes through them. 
This can be seen more clearly in figure an enlarge- 
ment of figure ([IJ . At each new H resonance caused by 
a downward fluctuation of the density the matter mixing 
angle ^13 will swing from a value close to it's high density 
limit of 90° through 45° towards its vacuum value and 
then back again: for an upward fluctuation that crosses 
the resonance density the matter mixing angle does the 
opposite swinging from close to the vacuum value up to- 
wards 90° and then down again. Either way, the non- 
adiabaticity parmeter (40[ - which describes the mix- 
ing between matter states 1^2 and - is proportional to 
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FIG. 2: An enlargement of figure U]) to show the detailed 
neutrino evolution in the MSW H resonance region. 



the derivative dOi^/dr^ . r23 is largest in the region of the 
resonance and if ever r23 ^ 1 then we have significant 
mixing between states 1^2 and and thus a change in the 
survival probability P22- Similar results are obtained for 
the antineutrino matter eigenstate survival probabilities 
Pii = P{j^i vi) and P33 = P{v3 —J- P3) in the case 
of inverted hierarchy because now it is the antineutrino 
non-adiabaticity ri3 which is proportional to the deriva- 
tive ddi^/dr and the H resonance mixes states Di and 

Figure ^ also allows us to observe the perturbative 
effect from non-resonant fluctuations. At non-resonant 
fluctuations the matter mixing angle 613 now varies over a 
smaller extent. The closer the fluctuation approaches the 
resonance the larger the variation. Despite the smaller 
variation of 613 at a non-resonant fluctuation compared 
to a resonant fluctuation the evolution of the neutrino 
state is governed by the derivative of this angle which 
can be large. For this reason we understand why non- 
resonant fluctuations can lead to the evolution of P22 as 
observed. 

The evolution of a neutrino with a particular energy 
and some set of oscillation parameters through the tur- 
bulent density profile is unique for each instantiation of 



^ More exactly one observes that both the and the Fia 

non-adiabaticity parameters defined in |40| are proportional to 
d6i3 / dr but as a consequence of the evolution of the matter mix- 
ing angle 612 in a normal hierarchy the non-adiabaticity function 
Fia is suppressed. 



5 



20 
15 

5 




0.5 



^ 0.5 



I ' I 



I ' I ' I 



0.2 



0.4 0.6 O.i 



FIG. 3: A normalized frequency histogram of 1012 calcula- 
tions of Pii (top panel) P22 (middle panel) and P33 (bottom 
panel) for E = 25 MeV neutrinos. The hierarchy is normal, 
sin^(26'i3) = 4 X 10~* and = 0.01. 
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the turbulence. If we change the F then we find a com- 
pletely different final state. If we repeat this exercise 
many times then we can construct an ensemble of final 
states from which we can study the general net effect 
of the turbulence i.e. the state of the neutrino after it 
has passed through the entire turbulence region and ex- 
its the supernova. The results of such an exercise are 
shown in figure ([3]) which displays the survival probabil- 
ity distribution of the three neutrino matter eigenstates 
in the case of E' = 25 MeV neutrinos, a normal hierarchy, 
sin^(26'i3) = 4 X 10""^ and = 0.01. Since mixing at the 
H resonance is between states U2 and z/3 and also the am- 
plitude of the fluctuation is relatively small the survival 
probability Pn is a delta-function at Pn = 1. Larger 
fluctuation amplitudes lead to broken HL factorization 
[33| . In contrast one can see that the distributions for 
P22 and P33 are consistent with uniform and the average 
of the ensemble for these two probabilities is 1/2. This 
implies that full 'depolarization' occurs for matter states 
1^2 and 1/3: one starts with a pure neutrino state and ends 
up with a flavor state that is completely random so that 
all information of the initial state is lost in the averages. 

As we decrease the fluctuation amplitude the number 
of new resonances created tends to zero and also the 
perturbative effect from non-resonant fluctuations disap- 
pears. But surprisingly we find that the turbulence docs 
not disappear quickly with C*. Turbulence continues to 



FIG. 4: A very fine energy resolution of Pn between 40 and 
41 MeV. The bottom panel is for no fluctuations, the next 
is C* = 10~^ then C^, = 10"^ and C* = 10"\ The results 
correspond to inverted hierarchy and a 'large' third neutrino 
mixing angle given by sin^(2^i3) = 4 x lO"**. 



have a large impact upon the final state probabilities one 
finds in an ensemble even for fluctuation amplitudes as 
small as ^ 10^^. We have investigated the reason 
for this extreme sensitivity and have found that it is due 
to a change in the relative phase difference between the 
extant resonances of the profile. That the sensitivity to 
small C* is due to the combination of phase effects and 
turbulence can be deduced from figure (|4]). Here the hi- 
erarchy is inverted but this makes no difference to the 
argument. Going from the lower to upper panels one can 
see that for values of C^, = lO^'^ — 10~^ the turbulence 
does not change the correlation energy scale of ^ 50 keV. 
Only for C^f > 10~^ does one observe a decrease in the 
correlation energy scale to ~ 10 keV as the turbulence 
effects start dominating over the phase effects. From the 
similarity of the middle panels of the figure with the case 
of no fiuctuations one concludes that the number of res- 
onances must be the same for all of them. Once this is 
appreciated one can then understand the sensitivity to 
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from the case of two resonance phase effects discussed in 
Kneller & McLaughhn [l^ and Dasgupta & Dighe [l6j . 
Using a two-flavor approximation of the mixing between 
states a and (3 the crossing probabihty for a neutrino 
passing through two resonances is 

Pc = Pi {I - P2) + (l - Pi) P2 

+2^PiP2{l-Pi){l-P2) cos$ (5) 

where Pi{ri) are the crossing probabihtics for the reso- 
nances separately and 

$ = / dr{K~kp) (6) 

is an interference term in which the r's are the locations 
of the resonances and the fc's are the matter eigenval- 
ues. Turbulence leads to fluctuations of ka — kp, rel- 
ative to the case without turbulence, and it also shifts 
the resonance positions ri and r2 ■ But if the fluctuation 
amplitude is small the changes in the resonance posi- 
tions and the gradient there are also small and so Pi and 
P2 are essentially unaltered. This leaves just the fluc- 
tuations of Phi due to the small shifts in the resonance 
positions and the difference of the eigenvalues ka — kp be- 
tween them. In order to generate a large effect we need 
a change 5^ ^ 1 but this is typically small compared to 
$ itself which is often $ 3> 1 because the resonances are 
very far apart. If (5$ ^ 1 then Pc will vary over the 
range AFc = 4^Pi P2 (1 - Pi){l - P2) and one may 
show that APc < 1/2 i.e. less than the spread of the 
uniform distribution shown above. 

So with this understanding of how neutrinos can be 
sensitive to such small amplitude fluctuations we present 
the full evolution of the survival and crossing probabili- 
ties as a function of for this case of a normal hierar- 
chy and 'large' Oi^. Figures ([5]) and ([6|) show plots of the 
probabilities Pn, P12, P22 and P23 as a function of the 
matter density fluctuation amplitude for 60 MeV neutri- 
nos in a normal hierarchy with sin^ 26113 = 4 x 10~*. The 
energy has changed compared to previous figures simply 
for the sake of clarity of the figure, the results are other- 
wise typical. And to demonstrate that the results are in- 
sensitive to the number of k modes we present results for 
Nk = 100 in figure © and in figure © Nk = 1000. From 
these figures four different ranges of can be identified. 
For values smaller than ^ 10~^° turbulence effects are 
negligible and the probabilities one derives are always the 
same and determined by the MSW resonances of the un- 
derlying adopted profile. For values between ^ 10~^" 
and 10~^ the turbulence effects and the phase effects in- 
duced by the multiple MSW resonances of the underlying 
profile, that must occur even in the absence of fluctua- 
tions, combine to produce broad variations over the en- 
semble but smaller than the depolarization limit. In the 
figures one observes that the rms variance of P22 and P23 
plateau at cr = 0.23, whereas the two flavor depolarized 
limit would give a — 0.28. In the region 10~^-10~^ tur- 
bulence effects start dominating over phase effects and 




FIG. 5: The mean and the variance of the neutrino matter sur- 
vival probabilities for a i? = 60 MeV neutrino, as a function 
of the noise amplitude. At every value of C* these quantities 
are calculated from at least 1000 samples of the turbulence F. 
The mean and rms variance of Pn is represented by the solid 
line in both panels, P22 is shown by the long-dashed, P12 and 
P23 are the short-dashed and dash-dot lines respectively. The 
hierarchy is normal and sin^(2Si3) = 4 x 10~^. The number 
of k modes is Nk ~ 100. 



lead to depolarization of states 1^2 and 1^3 but not vi. 
Note that figure ^ shows the distribution for the sur- 
vival probabilities when is in this range and one ob- 
serves in figures (O and ^ that a reaches 0.28 for these 
values of C+. Finally, for ^ 0.1 we start to transit 
to the case of three flavor depolarization where all three 
matter states start to mix due to broken HL factoriza- 
tion [3^. The distributions start to become triangular 
and the mean value drops towards 1/3. This transition 
to three flavor mixing is studied more carefully in the 
next section. 

Investigating how turbulence effects are modified when 
one decreases ^13 is important because the third neu- 
trino mixing angle may not turn out to be close to the 
present Chooz limit and able to be measured in the 
forthcoming reactor and first superbeam experiments. 
In a core-collapse supernova without turbulence, the 
non-adiabaticity of the H resonances of profiles from 
spherically symmetric simulations becomes large for val- 
ues of 013 smaller than the Dighe-Smirnov threshold of 
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FIG. 6: The same as figure © but with Nk = 1000. 
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FIG. 7: The survival probability Pu for a E = 25 MeV an- 
tineutrino as functions of the radius r but now with the value 
of ^13 set to sin^(2^i3) = 4 x 10^^. The instantiation of the 
density is shown in figure ([T]). Again, the turbulence region 
is identified by the vertical dashed lines and the evolution 
through the base profile as the dashed line. The hierarchy is 
inverted and the fiuctuation amplitude is C* =0.1. 



sin 2013 ^ 10~^ [42]. But when wc include turbulence 
the limiting behavior that survival probabilities for the 
states that mix at the H resonance should tend to zero 
as — >■ need not occur because the smallness of any 
non-zero can always be compensated by increasing 
the amplitude and/or extent of the matter fluctuations. 
However the fluctuation amplitude cannot be increased 
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FIG. 8; A normalized frequency histogram of 1687 calcula- 
tions of P22 for a i? = 25 MeV neutrino with a fixed under- 
lying profile and C* — 0.1. Here ^13 is set to sin^(2^i3) = 
4 x 10~^ and the hierarchy normal. 



beyond reason and the region where turbulence occurs is 
limited to the space between the shocks in our model. It 
is due to these 'upper limits' on the compensating vari- 
ables that indeed we can expect to reach the limiting 
behavior of non-adiabatic evolution for the H resonant 
mixing states as ^13 — )■ at some point but note it will 
not necessarily be at the same value as the Dighe-Smirnov 
threshold. This expectation is borne out when one inves- 
tigates the effect of turbulence for small values of 6113 
since one observes there can be significant, i.e. greater 
than 0(1%), changes in the neutrino survival probability 
for values of sin^ 26*13 = 4 x 10~^, as shown in figure ([7|). 
For this particular figure the hierarchy is inverted so the 
H resonance mixes the antineutrino states Di and for 
a normal hierarchy mixing is between i>2 and but the 
essential behavior is the same for both. 

As for the case of large ^13, each instantiation of F 
leads to a different final state and again we can con- 
struct the distribution of the possible final states by 
Monte Carlo methods. The distribution of the final state 
P22 for the case of a normal hierarchy, E ~ 25 MeV, 
sin^(26'i3) = 4 X 10"^ and = 0.1 is shown in figure 
([8|). We see that the distribution is not a delta- function 
at zero corresponding to non-adiabatic propagation but 
rather appears to have an exponential behavior. 

With figure ^ we summarize the entire behavior of 
the turbulence effect as a function of ^13 for a fixed mat- 
ter density fluctuation amplitude C* = 0.1. For 'large' 
angles - above the Dighe-Smirnov threshold- we are in 
the depolarized limit for this value of = 0.1, below the 
threshold we are in a scaling regime where P22 oc . Set- 
ting a limit that observable changes require {P22) > 0.01 
for present neutrino detectors then we see that turbulence 
induces observable sensitivity to ^13 for values above 
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FIG. 9: Mean value of the second matter eigenstate for 
E — 25 MeV neutrinos as a function of third neutrino mix- 
ing angle, in degrees, for the case of normal hierarchy and 
C* = 0.1. 



non-adiabaticity parameters and T13 which depend 
upon the derivatives of the matter mixing angles 9i2 and 
013 respectively. For both, the matter mixing angle varies 
between the high density limit of zero to its value in vac- 
uum at resonances. For 612 this occurs at L resonances, 
for 9i3 at H resonances. Thus we expect that the larger 
the value of the vacuum angle the greater the deriva- 
tive and the non-adiabaticity parameters T12 and T13 
and, therefore, the greater the mixing between the mat- 
ter eigenstates Di and D2 or Di and D3. The difference 
between the mixing of neutrinos states and vi-D^ 

is that the vacuum value of 612 is well known and large 
while the value of ^13 is neither large nor known. This 
effect of of non-resonant mixing between states Di and i>3 
is the same as found in Ref.[3l| where it is also shown 
that the turbulence impact on the non-resonant neutrino 
channel disappears for a decreasing value of ^13 though 
we would like to emphasize that our findings are for a 



0.02° or, equivalently, sin^ (26*13) > 4 x 10" 



B. Three flavor resonant and non- resonant mixing 

To complete our study of turbulence effects we look 
at large amplitude fluctuations, their breaking of HL 
factorization and the impact in the channels where no 
MSW resonance(s) occur. We take the case of normal 
hierarchy as an example so the non-L-resonant and non- 
H-resonant transitions occur for antineutrinos. In fig- 
ure (|I0|1 we present the turbulence density profile and 
the corresponding neutrino and anti-neutrino probabili- 
ties for the case of sin^ (26*13) = 4 x 10"'' and = 0.34. 
The amplitude of the fluctuations here is so large that 
HL factorization is broken. This can be seen from the 
evolution of Pn in the middle panel and noting that it 
occurs at the point where Ve in the top panel approaches 
the L resonance at r ^ 60,000 km. Also note that that 
breaking HL factorization does not depend upon the hi- 
erarchy. In the lower panel we see that simultaneously 
there is a sudden change in the anti-neutrino transition 
probabilities Pn and P12 at the point where the L res- 
onance occurs but we do not observe any variation of 
P13 which remains at zero for all r in this case. In or- 
der to observe non-trivial evolution of P13 we need to 
increase ^13 close to the Chooz limit. An example with 
much larger ^13 is shown in figure for the case of 
^'i3 = 9° and — 0.1. Again HL is broken, this time 
at r ~ 50, 000 km, causing a sudden decrease in Pn and 
simultaneous jumps in Pn and P12 but we also notice 
from the lower panel that P13 is already non-zero by this 
point; mixing between states Di and P3 had commenced 
when Ve approached the H resonance at r ^ 40, 000 km 
as shown in the top panel. Both these non resonant mix- 
ing effects can be understood through the antineutrino 
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FIG. 10: Impact of turbulence on tire first matter eigenstates 
for the neutrino survival probability (middle) and for the anti- 
neutrino survival and appearance probabilities (bottom). The 
results correspond to normal hierarchy, sin^(2^i3) = 4 x 10~*, 
= 0.34 and E = 25 MeV. The top panel shows the mat- 
ter density profile with noise and again the vertical dashed 
lines indicate the turbulence region while the upper horizon- 
tal line is the H resonance and the lower the L resonance for 
this particular neutrino energy computed using the two-flavor 
formula. 
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FIG. 11: Impact of turbulence on the first matter eigenstates 
for the neutrino survival probability (middle) and for the anti- 
neutrino survival and appearance probabilities (bottom). The 
results correspond to normal hierarchy, 613 =9°, =0.1 
and E = 25 MeV. The top panel shows the matter density 
profile with noise and, like figure pO|) . the vertical dashed 
lines indicate the turbulence region while the upper horizontal 
line is the H resonance and the lower the L resonance for 
this particular neutrino energy computed using the two-flavor 
formula. 



very different set-up than the one in |3l| . It is impor- 
tant to note that the non-resonant effect from 612 does 
not depend upon the hierarchy, it will always occur for 
the anti-neutrinos and never for the neutrinos because 
the sign of 6m^2 is known. In contrast the non-resonant 
effect from ^13 will switch to the neutrinos if we consider 
an inverted hierarchy. 

Thus we find two effects from turbulence in the non 
resonant channels. One from broken HL factorization 
causing non-L-resonant transitions and the other from 
the constant transitions from the high density to vac- 
uum values of 6*13 . As with the case of small 6*13 and the 
H resonant channel, the distributions of the non resonant 
transition probabilities for an ensemble have an exponen- 
tial like behavior. An example of the distributions for 
the particular case of for a normal hierarchy, ^13 = 9°, 

= 0.2 and E = 25 MeV is shown in figure P^ . 
As explained above, the average values of the distribu- 
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FIG. 13: Variation of the mean anti-neutrino survival prob- 
abilities (-P11) (solid), (-P12) (short dashed) and (-P13) (long 
dashed) as a function of the vacuum mixing angle ^13 for the 
case of a normal hierarchy, C* = 0.1 and E — 25 MeV. The 
symbols indicate the values of ^13 where averages were calcu- 
lated. 



tions scale differently with ^13 and the variation of the 
turbulence effect with the third neutrino mixing angle is 
shown in figure p^ . While the first matter eigenstate 
survival probability and appearance probability to the 
second matter eigenstate (P12) are constant the mixing 
with the third matter eigenstate, {P13), scales as 6*^3 and 
so disappears as 6*13 decreases. We also observe that the 
mixing between states j>i and j/3 surpasses that between 
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FIG. 14: Variation of the mean anti-neutrino survival prob- 
abilities (-Pii) (solid), {P12) (short dashed) and (-P13) (long 
dashed) as a function of the fluctuation amplitude C* for a 
vacuum mixing angle 6^13 = 9°, the case of a normal hierar- 
chy, and E — 25 MeV. The symbols indicate the values of C* 
where averages were calculated. 



Pi and P2 only for ^13 > 3°. Finally, the behavior of 
(Pii), (P12) and (P13) as a function of the fluctuation 
amplitude is shown in Fig. p^ . At small amplitudes 
both (P12) and (As) scale as but as increases (P12) 
appears to increase its sensitivity to this parameter. 



IV. CONCLUSIONS 

We have performed a detailed investigation of turbu- 
lence effects after solving the evolution equation for neu- 
trinos passing through turbulent density profiles created 
by adding matter density fluctuations to a ID density 
profile obtained from hydrodynamical simulations. Our 
treatment differs from most previous studies since we do 
not solve the neutrino evolution equations for the aver- 
aged probabilities but rather for the neutrino amplitudes. 
The results have highlighted new turbulence effects for 
that period when the shocks are in the region of the H 
resonance for neutrino energies of order 0(10 MeV) and 
can be summarized as follows. 

• For large third neutrino mixing angle but fluctu- 
ation amplitudes below ^ 0.1 we find turbu- 
lence effects only appear in the channel with the 
H resonance i.e. the neutrinos for a normal hi- 
erarchy, the antincutrinos for an inverted hierar- 
chy. More specifically wc find that if the matter 
density fluctuation amplitudes are tiny, no turbu- 
lence effects appear as expected; for intermediate 
values, 10^^° ^ ^ 10~^ the neutrino matter 
survival probabilities are dominated by the phase 
effects between the pre-existing resonances of the 
profile but with fluctuating relative phase that are 
sufficient to produce a wide variation of the final 



states; for larger values of the matter density fluc- 
tuation amplitude, 10~^ ^ ^ 10~^, turbulence 
effects begin to dominate and lead to two-flavor de- 
polarization. The sensitivity to amplitudes as small 
0.001% would apparently indicate that turbulence 
effects in this channel arc inevitable but, at the 
same time, obscrvationally turbulence effects will 
be difficult to distinguish from phase effects once 
one takes into account time and energy binning of 
a signal, as pointed out in Calais et al. f2l| . 

• When the matter density fluctuation amplitude is 
large, > 0.1, we have found a sensitivity in the 
H resonant channel to the third neutrino mixing 
angle below the Dighe-Smirnov threshold down to 
sin2(26li3) > 4 X 10"^. 

• Large amplitude matter density fluctuations, ^ 
0.1, lead to broken HL independent of ^13 and the 
hierarchy. This leads to resonant mixing between 
states vi and 1^2 cmd non-resonant mixing between 
states i>i and i>2. 

• Once again, for large amplitude matter density fluc- 
tuations, ^ 0.1, we find that non- resonant mix- 
ing may occur due to the fiuctuations of the matter 
mixing angle ^13. The degree of mixing is a func- 
tion of 6*13 with a mean value that scales as Of^. The 
channel where this effect appears depends upon the 
hierarchy appearing in the antincutrinos for a nor- 
mal hierarchy where it leads to mixing between Pi 
and P3 and in the neutrinos for an inverted hierar- 
chy where it causes mixing between states 1^2 and 

When the matter density fluctuations are small turbu- 
lence can only cause two flavor mixing. When the ampli- 
tude is large three flavor mixing may occur in the neu- 
trinos and antincutrinos simultaneously but, depending 
upon the hierarchy and ^13, not necessarily to the same 
degree. Let us summarize our flndings in this case of 
large amplitude fluctuations: 

• In the case of normal hierarchy: for ^13 above the 
Dighe-Smirnov threshold strong mixing occurs in 
the neutrinos and 3-flavor depolarization may be 
approached. On the other hand, if ^13 is close to the 
Chooz limit, 3-flavor mixing of the antineutrinos 
occurs, due to the two non-resonant effects; oth- 
erwise only two-flavor mixing takes place because 
of the non-resonant effect from 612 fluctuations. 
When 9i3 is below the Dighe-Smirnov threshold, 
the neutrinos still exhibit three flavor mixing due 
to broken HL factorization and the residual sensi- 
tivity to 013, that scales as Of^. 

• In the case of inverted hierarchy: for large fluctu- 
ations and 9i3 close to the Chooz limit, 3-flavor 
mixing occurs in the neutrinos, due to broken 
HL factorization, and the non-resonant effect from 
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6*13 fluctuations which will turn into 2-flavor non- 
resonant mixing of vi and 1^2 as ^13 decreases. 
Three flavor mixing of the antineutrinos takes place 
for 013, above the Dighe-Smirnov threshold, but it 
is unlikely that depolarization is reached because 
the mixing between Di and D2 is non-resonant. For 
9i3 below the Dighe-Smirnov threshold, some resid- 
ual three-flavor mixing will remain, but it will be 
small, because now the two turbulence effects are 
perturbative in this channel. 

Several aspects related to turbulence remain open. 
While our treatment of turbulence can already be con- 
sidered, in many respects, realistic it is clear that fu- 
ture studies are awaited where matter density fluctua- 
tions are drawn from successful core-collapse supernova 
simulations at many different epochs. Also, until we can 
calibrate these one-dimensional studies the actual turbu- 
lence effects we should expect is unknown. 

Lastly, one important question that needs further in- 
vestigation is the extent of the 'smearing' or 'blurring' of 
turbulence features in a burst signal. Smearing can oc- 
cur due several reasons. First there is the smearing due 
to the finite size of the source. The neutrinos we shall 
receive from the next Galactic supernova will not travel 



along a single line of sight so the neutrino flux at a given 
arrival time and neutrino energy (ignoring time of flight 
issues [i^l) is 'averaged' over a column of ~ 20 km i.e. the 
apparent diameter of the proto-neutron star source. Neu- 
trinos emitted from different locations upon the source 
will pass through different small scale, < 20 km, fluctua- 
tions but the turbulence on large scales will be the same 
for all the neutrinos we receive in the burst. Second, 
a smearing in energy will occur in the detector due to 
the intrinsic energy resolution of the detection reaction 
and the extrinsic resolution of the detector. If the com- 
bined effect of intrinsic and extrinsic resolutions leads to 
a smearing over an energy range a then features on scales 
smaller than cr will be unobservable. Finally, smearing 
in both time and energy will also be introduced because 
of the necessity of binning the detected events. As with 
the detector response, short time scale changes in the 
flux or features with small energy widths will be smeared 
reducing their observability. Studying how all these ef- 
fects impact the observability of the turbulence effects 
pointed out in the present paper is not trivial and will 
be the object of future investigations. [4lj |. 
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